Full counting statistics of spin transfer through ultrasmall quantum dots 
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We analyze the spin-resolved full counting statistics of electron transfer through an ultrasmall 
quantum dot coupled to metallic electrodes. Modelling the setup by the Anderson Hamiltonian, we 
explicitly take into account the onsite Coulomb repulsion U. We calculate the cumulant generating 
function for the probability to transfer a certain number of electrons with a preselected spin orien- 
tation during a fixed time interval. With the cumulant generating function at hand we are then 
able to calculate the spin current correlations which are of outmost importance in the emerging 
field of spintronics. We confirm the existing results for the charge statistics and report the discov- 
ery of the new type of correlation between the spin-up and -down polarized electrons flows, which 
has a potential to become a powerful new instrument for the investigation of the Kondo effect in 
nanostructures. 
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Modern microelectronics is one of the most success- 
ful technologies ever conceived by humankind. However, 
in recent years a number of limitations, which can slow 
down or even stop further progress began to come to 
the fore. One possible way to overcome these difficulties 
is to switch from charge current processing to spin cur- 
rent and spin configuration processing. Their advantages 
are so enormous that recently a completely new scientific 
field of spintronics has been established [l|, l2| . 

In the conventional microelectronics the properties of 
the basic circuitry elements are characterized by a num- 
ber of different quantities - by the nonlinear current- 
voltage relations, by the current noise spectra, by the 
current correlations of third order (third cumulant) etc 
0,0]. However, there is one characteristic which (at least 
in the low frequency range) contains information about 
correlations of all orders. This is the so-called full count- 
ing statistics (FCS), which answers all above questions 
by providing the probability distribution P(Q) to trans- 
fer a certain amount of charge Q during a waiting time 
interval T 0,0]. 

While there is by now a vast amount of literature avail- 
able on the charge transfer statistics, its spin-resolved 
relative remains to a larger part unknown with some no- 
table exceptions 0, 0, @, E3]. We would like to close 
this gap and analyze the combined statistics of spin and 
charge transfer through ultrasmall quantum dots with 
genuine repulsive interactions. The distinctive feature 
of these devices are their extremely small lateral dimen- 
sions which allow for only very few energy levels to take 
part in the transport processes. Typical realizations are 
nanoscale heterostructures on the semiconductor basis, 
or individual molecules coupled to metallic electrodes, 
which are even smaller [TJ 1 12L IT3L [3] • Both types of 
systems became available only during the last decades 
and are promising candidates to become basic building 



blocks of future nanoelectronics and spintronics circuitry. 

The archetype system to describe such structures is 
the single-impurity Anderson model [HEl. It consists 
of a local fermionic level (also called dot level) which 
is filled or emptied by d\ , d„ creation and annihilation 
operators and N different electronic continua modeling 
the electrodes via fermionic fields ipao-ix), 

N 

a— ± a— 1 cr=± 

where A+ah =: b a is responsible for the energy of the dot 
in the magnetic field h = §^5/2, when it is occupied by 
an electron with spin orientation a. The dot level energy 
A is an additional parameter which can in practice be 
changed by varying the voltage on the background gate 
electrode. The coupling of dot and leads is achieved by 
a local tunneling contribution 



H T = Y^^ a [4^w(z = 0) + h.c] 



(2) 



with energy independent tunneling amplitudes j a . In 
addition to these terms one has to take into account the 
electrostatic repulsion 



H v = U(d\df - n )(d\di - n ) , 



(3) 



reflecting the energetic cost U of double dot occupation 
with respect to the symmetric value no = 1/2. Due to 
this normalization, the particle-hole symmetric case cor- 
responds to A = 0. The full system Hamiltonian is the 
sum of all three contributions H = Hq + Ht + Hjj. 

The technology for the calculation of the FCS is by 
now far advanced and allows for a number of different 
approaches. In the most widespread one, the quantity 
of interest is the so-called cumulant generating function 
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(CGF) ln X (A i; A T 



lnx(A 



17, m 



Its succes- 



sive differentiation with respect to the counting fields \ a 
yields the respective irreducible momenta {{5Q™}} for the 
probability distribution to transfer SQ a charges with spin 
orientation a through the system during the waiting time 

r, 



((5Q n J = (-i) n ^ln X (Xa)\ x ._ 



(4) 



In analogy to the approaches taken in ljj 20| . the first 
step in the calculation of the CGF is to endow the tun- 
neling Hamiltonian ((2|) with counting fields X a . Outside 
of the waiting time interval < t < T the counting fields 
are zero. The tunneling Hamiltonian then transforms to 



Hn 



[e* W2 4Vw+h.c. 



(5) 



In the noninteracting case (U = 0, resonant level model 
20l . 2lj ) this quantity can easily be calculated by resum- 
mation of the perturbation series in j a or by applying 
the Levitov-Lesovik formula [5|. Assuming the measure- 
ment time T to be large such that switching effects can 
be neglected, one finds (setting e = h = 1) 

ln Xo (A CT ) =TJ2[ ^{l + Y.Tf^Wil-np) 

a J 71 a/3 

x [ e ^--^)-l]|, (6) 

where np(u>) denotes the Fermi distribution in lead (3. 
The energy dependent transmission coefficients are given 
by 
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a 1 /3 



(^-(5 CT ) 2 + r 2 



(7) 



where = irwyp (with the density of states v at the 
Fermi level) is the hybridization of the dot with lead (3 

andr = ^ /3 r /3 . 

While the non-interacting result ([6]) was derived for an 
arbitrary number N of fermionic leads, we shall restrict 
ourselves henceforth to the symmetric two level case, i. e., 
we assume N = 2 and Tl = T^. The chemical potentials 
of the leads are assumed to be at hl,r where V = (J-l—I^r 
denotes the applied voltage. 

It is quite inefficient to calculate x(\ a ) using an ad- 
ditional expansion in U in the same way. As has been 
realized in Ref. [22| in a different context, as long as one 
is only interested in small U, the CGF can be calculated 
by a simple linked cluster like calculation. Thereby the 



full A CT dependence can be shifted onto the unperturbed 
Keldysh Green's functions D(u>) of the dot level. Using 
the notation of [23j , these are given by 

D (lo) = [(w-^)+532*r /J (n^-l/2)]/X>(w), 
D+- («) = 2iY p e~ lX ^ ( nfl - 1)] /X»(w) , 



= -[D—(u)]', 
where we defined 

V{uj) = (oj~8 a ) 2 + T 2 

+ 4y^r a r /3 n a (l - np) 

a/3 



(8) 



(9) 



Note that the presence of the counting field causes a vi- 
olation of Keldysh's sum rule. Using a linked cluster 
expansion, the exact CGF \nx(X a ) can be expressed as a 
correction to the noninteracting CGF ([6]). The quantity 
we have to evaluate is then 



x(A<x) = xo(A CT )(T c exp [■ 



i / dtHu(t)]). 
Ic 



(10) 



The expectation value is to be taken with respect to the 
noninteracting ground state and therefore contains the 
A-dependent Green's functions (JSJ. It may appear that 
since the lowest order expansion in U only contains the 
dot occupation numbers n a , one is not expecting any 
counting field dependence to survive. However, this is 
no longer valid in the case of explicitly (quite artificially) 
time-dependent A CT . In the limit 



max{A, h, V}/T < 1 , 
the first order contribution is given by 



mx (i)(A CT ) = -^J^^' 



e~ lX " - 1) 



(11) 



(12) 



where a = —a. As it depends on 8 a , this term only con- 
tributes for finite magnetic field (h =/= 0) and/or broken 
electron- hole symmetry (A ^ 0). 

The second order contribution is given by two different 
diagrams, see Fig. [TJ One of these is again proportional 
to the magnetic field and contains the average dot oc- 
cupation numbers while the other one is the double shell 
diagram. The calculation procedure is rather lengthy but 
straightforward and results in the following CGF expan- 
sion m the limit (fTTj) , 
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FIG. 1: One first order (a) and two second order (b,c) dia- 
grams. 



where we introduced the equilibrium even/odd suscepti- 
bilities (correlations of nj with nj and nj,, respectively), 
which are known to possess the following expansions for 
small U 



Following the reasoning along the lines of Ref. [2fJ, we 
may speculate that (fT3|) is the exact result to all orders 
of U and T as soon as one inserts the exact values for 
Xe,o which have been obtained by, e.g., Bethe ansatz cal- 



culations [27|, l28( . Now we are in a position to establish 
contact to known results and to discuss new effects. Thus 
far, similar results have been obtained only for the charge 
transport statistics of the same system for the much more 
restrictive particle-hole symmetric parameter constella- 
tion [2(J. The complete spin resolved statistics for large 
transmission T is given by the following CGF, 
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TV 

Inx(Xa) = iG TV(\i+\ l ) + — -J2( X eS, + Xo6„) 2 (e-^ - 1) 
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where Go = l/(2n) is the conductance quantum per spin 
orientation. In order to go over to the pure charge CGF 
one has to set Aj = Aj = A. In this case, the result 
of Ref. [2(| is perfectly reproduced for h = A = 0. 
Moreover, it had been speculated that while the terms 
containing a single A correspond to single electron tun- 
neling events, the term with the doubled counting field 
is brought about by a coherent electron pair tunneling 
in a singlet state. Eq. ([T5)) represents the proof of this 
conjecture since the term giving rise to 2A part indeed 
stems from the contribution originally containing the sum 

A T + A X . 

Yet another justification of the validity of (| 1 5[) for arbi- 
trary U is brought about by comparing the above result 
to the spin-resolved statistics of charge transfer through 
a Kondo impurity in the unitary limit presented in [221 ] . 
Similar to the parameter mapping for the conventional 
current statistics one identifies the two limits of small 



and large U (Kondo regime) by 

HTk=Xo/T and a/T K = Xe /T, (16) 

where Tk is the Kondo temperature and (f> and a are 
Fermi liquid parameters of the Kondo fixed point [2^ |. 
The similarity of these two results can be traced back to 
the similarity of the corresponding Hamilton operators, 
which not only both contain a resonant level part but 
also possess analogous interaction terms. 

Next, we would like to discuss the linear response (lin- 
ear in V) contribution. It can easily be verified that un- 
der the conditions (jTTJ) one obtains the following CGF, 

lnx(A .)h:„ = (17) 
This result perfectly coincides with the conjecture of the 
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binomial theorem formulated in [20[ . It predicts that the 
linear response charge transfer statistics of any interact- 
ing region coupled to noninteracting continua is binomial 
and governed by the value of the transmission coefficient 
at the Fermi edge. In fact, the road to the construction 
of the spin-resolved CGF from Eq. (26) of [U (which 
contains only the charge transfer generating function) is 
very natural and intuitive: the logarithms with different 
signs in front of the magnetic field term should contain 
counting fields for different spin projections. 

The spin current statistics can easily be recovered from 
the above results after transition to the charge current 
and spin current counting fields A, // via Afj = A ± /z. 
One feature of (fT5|) is the fact that the odd cumulants of 
spin currents are only non-zero in finite field and for the 
particle- hole asymmetric case A ^ 0. This is the precise 
condition for the spin flow existence in a conventional 
noninteracting resonant level system as well. In the linear 
response regime the corresponding odd order cumulants 
are then given by 

«T - S Ql ) 2n+1 )) = - 2J ^ (Xl ~ xl) • (18) 

The even order cumulants arc non-universal but n- 
independent as well, so that the ratio of even/odd orders 
(it can be seen as a generalization of the Fano factor) is 
given by 

«Wt - SQ L ) 2n )) 



<(w T - 

A 2 (Xe 



\2n+l\ 



+ Xof + h 2 { X e-Xof 



(19) 



2Afc( X 2 - xl) 

Going beyond the linear response regime, we find that 
the most fundamental feature emerging from (|15[) is the 
existence of the invariant cross-cumulant 



((SQ\ l SQf)) = (-i) 



n-f m 



= (-1) 



7i-\-7n 
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lnx(A CT ) 



(20) 



for n, m > 1. Not only is this quantity non-zero in inter- 
acting systems only, it is also independent of magnetic 
field strength and (up to the sign) of its orders n, m. It 
exists in the strong coupling Kondo case as well and is 
found using (|16[) for the parameter translation between 
weak and strong coupling. Despite the formally identical 
mathematical shapes, the origin of this phenomenon is 
completely different for small U and in the Kondo regime. 
While in the weak coupling case it signifies the begin- 
ning of the spin singlet formation, in the strong coupling 
limit it starts to appear as soon as it becomes possible to 
break up (though virtually) the Kondo spin singlet. The 
amplitude of these correlations grows as one approaches 
the strong coupling fixed point. In principle, in addi- 
tion to the conventional linear conductance (which ap- 
proaches the unitary limit of almost perfect conductance) 



the cross-cumulant can also be regarded as a measure of 
how deep in the Kondo regime the system in question is 
being. 

To conclude, we have analyzed the non-equilibrium 
spin resolved FCS of the Anderson impurity model by 
calculating the CGF of the probability distribution to 
transfer a fixed amount of charge with preselected spin 
orientation during very long waiting time interval. Our 
results perfectly agree with existing predictions for the 
statistics of charge transfer. The CGF indeed supports 
the interpretation that the electron transport in such a 
system is mediated not only by single charge tunneling 
but by correlated transfer of electron pairs in a singlet 
state. Moreover, the emerging expressions confirm the 
previously conjectured statistics in finite field beyond the 
particle-hole symmetric situation. In the linear response 
regime, it turns out to be binomial and to factorize in dif- 
ferent spin channels. Finally, we have discovered a new 
type of correlation between the spin-up and spin-down 
currents: a cross-cumulant. It is universal and field in- 
dependent. In our view, it has a potential to become one 
of the quantities to measure and control the 'quality' of 
the Kondo effect in nanostructures. 
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